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Abstract 

We construct a flavor symmetry model based on the tetrahedral group A/± in which the right-handed 
neutrinos from the seesaw mechanism can be both keV warm dark matter particles and eV-scale 
i-^ \ sterile neutrinos. This is achieved by giving the right-handed neutrinos appropriate charges under 

the same Froggatt-Nielsen symmetry responsible for the hierarchy of the charged lepton masses. 

P-i 

We discuss the effect of next-to-leading order corrections to deviate the zeroth order tri-bimaximal 



mixing. Those corrections have two sources: (i) higher order seesaw terms, which are important 
when the seesaw particles are eV-scale, and (ii) higher-dimensional effective operators suppressed 
by additional powers of the cut-off scale of the theory. Whereas the mixing angles of the active 
neutrinos typically receive corrections of the same order, the mixing of the sterile neutrinos with 
the active ones is rather stable as it is connected with a hierarchy of mass scales. We also modify 
an effective model to incorporate keV-scale sterile neutrinos. 
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I. INTRODUCTION 



Apart from the direct proof of physics beyond the Standard Model (SM) in the form 
of neutrino masses lj, a somewhat more indirect proof is the presence of Dark Matter 
(DM) 3|. One can take the point of view that these two aspects of physics beyond the SM 
are connected with each other, i.e. that neutrino mass and Dark Matter are linked. We will 
assume this connection in the present paper. 

The most direct such relationship would be realized if the light massive neutrinos whose 
oscillations we observe in the lab are the DM particles. However, they would be Hot Dark 
Matter, and cosmological data is compatible only with a very small component of this 
form of DM, which in fact allows one to set limits on neutrino mass jij]. Typically the 
DM is assumed to be of the Cold Dark Matter (CDM) type, for which a WIMP (weakly 
interacting massive particle), as predicted in many supersymmetric theories, is the most 
popular candidate. However, Warm Dark Matter (WDM) is another possibility compatible 
with observations, and in fact could solve some of the problems of the CDM paradigm, in 
particular by reducing the number of Dwarf satellite galaxies or smoothing the cusps in the 
DM halos. At this point one should note that a sterile neutrino with mass at the keV scale 
and with small mixing to the active nentrinos is a WDM candidate if a mechanism! | 
to generate the correct amount of relic population is present 1 . See the reviews |8l-ll0| for 
summaries of mechanisms and the status of keV sterile neutrinos as DM. 



Sterile neutrinos heavier than the active ones are an ingredient of the seesaw mecha- 



nism 
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151 ] , whose existence is strongly hinted at from the fact that active neutrino masses 
are extremely small. Here, however, the right-handed neutrinos are "naturally" of order 10 10 
to 10 15 GeV, and if one wishes to make one of them a WDM candidate one has to arrange 
for this mass to come down to the keV level. The following possibilities exist: 

• theories with extra dimensions can exponentially suppress fermion masses, by localiz- 
ing them on a distant brane, for instance. This has been proposed to generate seesaw 



neutrinos of ke V scale in [2] , see also [l7j ; 



flavor symmetries jl8l . can predict that one of the heavy neutrino masses is zero. 
Slightly breaking this symmetry generates a neutrino with much smaller mass than 



the other two, whose masses are allowed b y t 
generate seesaw neutrinos of keV scale in 120 



r e sy mmetry. This has been proposed to 



2lJ, see also |22j; 



while the commonly studied flavor models with non-abelian discrete symmetries cannot 
produce a non-trivial hierarchy between fermion masses, the Froggatt-Nielsen mecha- 
nism is capable of this 23j. This has been proposed to generate seesaw neutrinos of 



1 keV sterile neutrinos could also provide an explanation for pulsar kicks 
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keV scale in [24j, see also [25] 



extensions or variants of the canonical type I seesaw often contain additional mass 
scales, which can be arranged to generate keV-scale particles. This has been proposed 
to generate seesaw particles of low scale in [26[, see also [25]. 



Note that both the Froggatt-Nielsen and extra-dimensional approaches require that the 
three right-handed neutrinos cannot be identified as a triplet of a flavor symmetry, which 
is very often the case in flavor symmetry models (see for instance the classification table 



for A 4 models in Ref. 27[). Furthermore, in that case there is no overall effect on the 
leading order seesaw formula M D /Mi, with being the right-handed neutrino mass and 
Mr, the Dirac mass. Both mechanisms will suppress Mi quadratically, while Mo is linearly 
suppressed, and hence their combination M D /Mi is left constant. 



In this paper we will apply the Froggatt-Nielsen mechanism to bring one of the heavy 
neutrinos from its "natural" scale down to the keV level. We will construct an explicit 
flavor symmetry model based on the group A±. As in many such models, there is also a 
Froggatt-Nielsen symmetry £/(1)fn to generate the observed hierarchy of the charged lepton 
masses; we will use this very same £/(1)fn for creating a WDM candidate from the heavy 
neutrinos. 

In addition, it should be noted that when one goes from, say, 10 15 GeV = 10 24 eV 
down to keV = 10 3 eV, it is not a big problem to reduce the mass by another 3 orders of 
magnitude. In this way one has generated one (or more) sterile neutrino(s) of order eV. This 
would be very welcome to explain long-standing issues in particle physics, astrophysics and 
cosmology. Those are the apparent neutrino flavor transitions at LSND and MiniBooNE, 



which together with the "reactor anomaly" 28|, |29| point towards oscillations of eV-scale 



sterile neutrinos mixing with strength of order 0.1 with the active ones (see Refs. |30l. |31| for 
recent global fits 2 ). In addition, several hints mildly favoring extra radiation in the Universe 



34 
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have recently emerged from precision cosmology and Big Bang Nucleosynthesis 
This could be any relativistic degree of freedom or some other New Physics effect, but has a 
straightforward interpretation in terms of additional sterile neutrino species. Although some 
tension between the neutrino mass scales required by laboratory experiments and the Hot 
Dark Matter limits exists within the standard ACDM framework, moderate modifications 
could arrange for compatibility 38[. Finally, active to sterile oscillations have been proposed 
to increase the element yield in r-process nucleosynthesis in core collapse supernovae (which 
seems to be too low in standard calculations, see e.g. [39|, |40|). It is rather intriguing that 
indications of the presence of eV sterile neutrinos come from such fundamentally different 
probes. 



2 Sometimes the result of the calibration of Gallium solar neutrino experiments 32[ is interpreted as the 
"Gallium anomaly" and is considered to be an effect of sterile neutrinos 33 1. 
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We note that a particul ar p henomenological model, the i/MSM [y Minimal Standard 
Model), has been proposed 41], in which one of the seesaw neutrinos is keV and the other 
two can generate the baryon asymmetry of the Universe either via leptogenesis (if they are 
heavy) or via oscillations when they have masses below the weak scale and are degenerate 
enough jg), |4l|. The idea to exploit the neutrinos of the canonical seesaw mechanism to 
account for the required eV and keV particles has been discussed in Ref. [42|, for instance. 
Here we provide a reasoning for the low-lying scales and add to the framework a flavor 
symmetry that yields at leading order tri-bimaximal mixing (TBM). In addition we modify 
an existing effective model, which does not contain the seesaw mechanism, by adding 
a sterile neutrino. Again, applying appropriate Froggatt-Nielsen charges gives the correct 
charged fermion mass hierarchy and a WDM particle. 

As a starting point in our seesaw models, we will leave the Froggatt-Nielsen charges of the 
seesaw neutrinos free, except for the one which is doomed to be the keV WDM particle. By 
properly choosing the charges of the other two, we can make one or two to be of eV scale, or 
keep both heavy (below or above the weak scale). Different and testable phenomenology in 
terms of short-baseline oscillations or neutrinoless double beta decay (Oz//3/3) is then present 
and characteristic for each scenario. For instance, if all neutrinos are below the momentum 
scale 100 MeV of double beta decay, the effective mass on which the amplitude depends 
cancels exactly. This is in contrast to the usually considered analysis of sterile neutrinos in 
double beta decay j25| (or our effective A 4 model), in which sterile neutrinos are simply added 
to the three active ones and treated as independent entities. Interestingly, this cancellation 
happens pairwise in our particular model, because the columns of the Dirac mass matrix 
are proportional to the columns of the lepton mixing matrix and each of the right-handed 
neutrinos is responsible for generating one light active mass. If this right-handed neutrino is 
lighter than 100 MeV, then its contribution to double beta decay cancels exactly with that 
of the associated light active neutrino. 

We take particular care in evaluating next-to-leading order (NLO) corrections to 
the model, which lead to deviations from TBM. Two sources for those corrections are 
considered. If the right-handed neutrino mass is ~ eV instead of the natural value 10 10 



to 10 15 GeV, then NLO seesaw terms 43l445j can be important. This is because the 



seesaw formula goes like Mp/Mi (1 + M^/Mf). It is easy to see that if Mi ~ eV and if 
M 2 D I Mi ~ 0.1 eV, M D should be around 0.3 eV, and hence the NLO seesaw term M 2 D jMf 
can generate effects in the percent regime. Another, more commonly studied source of 
NLO corrections stems from higher- dimensional operators suppressed by additional powers 
of the cut-off scale of the theory. The relative magnitude of those terms also depends on 
details of the model and of the scales chosen for the neutrinos and other particles. We 
show in particular that values of U e s compatible with recent fits 46|, |47[ can be obtained in 
our models. An important aspect is that the three mixing angles of the active neutrinos 
typically receive corrections of the same order, as is generically the case in flavor models. 
However, the mixing of the sterile neutrinos with the active ones is rather stable as it is 
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defined as a hierarchy of mass scales, thus stabilizing, for instance, the small mixing of the 
WDM neutrino with the active ones. 



The remaining parts of this work are organized as follows: in Section [TT] we present some 
model-independent features of the seesaw mechanism and its resulting phenomenology in 
the case that one or more of the right-handed neutrinos is light. Section II III introduces 
a seesaw model based on the A4 flavor symmetry, in which one of the three right-handed 
neutrinos acts as the WDM candidate. Various cases for the mass scales of the other two 
neutrinos are discussed, phenomenological consequences are figured out in detail, and the 
role of higher-order corrections is studied. Details of the NLO terms are delegated to the 
Appendix. Section IIVI details an effective theory with a single keV sterile neutrino added 
to an existing A4 model. Higher-order corrections and possible deviations from the exact 
TBM pattern are also discussed. We summarize and conclude in Section [V] 



II. LIGHT STERILE NEUTRINOS IN TYPE I SEESAW 



Before describing a specific model, we address the role of light sterile neutrinos in the 
type I seesaw, in particular the effect of NLO seesaw corrections to neutrino mixing param- 
eters as well as phenomenological consequences of light sterile states. 



A. NLO seesaw corrections 



In the canonical type I seesaw mechanism, one extends the SM particle content with 
three right-handed neutrinos (vf, v^) together with a Majorana mass Mr. The full 6x6 
neutrino mass matrix in the basis (u e , u^, u T , v\ y z/|, v^) reads 

Mr = ( ° T , (1) 

[Mr Mr) 1 ' 

where Mp denotes the Dirac mass term, and we use the LR convention for the Lagrangian. 
Assuming Mr S> Md, this mass matrix can be approximately diagonalized using a 6 x 6 
unitary matrix as 

U „ (l-&Bt B \ f Vv ()\ 
v -\ 1 - \&B )\0 V R )' 1 ] 

where B = MpM^ 1 + O [M D (M R ~ ) 3 ) governs the effect of higher-order seesaw correc- 
tions 43|-|45|. The matrix V u is given by 

M v = -M d Mr X MI = V v diag(mi, m 2 , m 3 )Vj (3) 
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with rrii being the light neutrino masses, and Vr diagonalizes the right-handed neutrino 
mass matrix, i.e. M R = V R diag(M 1; M 2 , M 3 ) Vg. 

In the ordinary type I seesaw framework, Mr is commonly chosen to be close to the 
Grand Unification scale (e.g. M R ~ 10 14 GeV), while M D ~ 100 GeV, so that the light 
neutrino masses are suppressed at the sub-eV scale, i.e. O^M^/Mr) ~ 0.1 eV. Therefore, 
the NLO seesaw corrections governed by B can be safely neglected. In models with keV-scale 
sterile neutrinos these corrections are also negligible. However, for models with right-handed 
neutrinos located at very low-energy scales, i.e. at the eV scale, Mr, ~ 0.1 eV is required in 
order to generate light neutrino masses. In that limit the NLO seesaw terms are significant, 
and B ~ 0.1 may lead to sizable corrections to neutrino mixing parameters. In the remaining 
parts of this work we will keep the NLO seesaw corrections up to 0(B 2 ). Note that the 
block-diagonalization of Eq. ([1]) by Eq. ([2]) is still approximately valid, as the remaining 
off-diagonal terms BM^/Mr ~ M L ,jM\ ~ 0.001 eV are much smaller than the 0(1) eV 
mass difference between the active and sterile neutrinos 3 . 



B. Active-active and active-sterile mixing 

In the basis where the charged lepton mass matrix is diagonal, the light active neutrinos 
mix via the 3x3 matrix (1 — \BB^V V , whereas the mixing between the active neutrino u a 
(a = e, r) and the sterile neutrino uf (i — 1, 2, 3) is given by 

9ai = [U u ]a,3+i = [BV R ] 

where M^ 1 = diag(M 1 _1 , , M^ 1 ). This illustrates that active-sterile mixing is defined as 
the ratio of two scales, Md and Mr. The interaction between each sterile neutrino uf and 
the entire active sector is 

E i*-i 2 - ( 5 ) 

a=e,fi,r 

In the setup described above with eV-scale right-handed neutrinos and Mr, ~ 0.1 eV, 
9i = OiMrtjMi) ~ 0.1 is obtained, which could provide an explanation for the short-baseline 
anomalies. In the same way, for a keV-scale particle and the same Dirac scale Mjy — 0.1 eV, 
one gets 9i ~ 10 -4 (see the discussion in the following subsection). 

With the above notation, it is not difficult to check that the standard seesaw formula in 



In our numerical calculations we did not use any approximations, but rather numerically diagonalized the 
full neutrino mass matrix, obtaining results consistent with the analytical calculations. 



M D {V^M^)V R 



[1V1 D V 1 



Rloa 



Ma 



(4) 
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Eq. ([3]) can be re-expressed as 



[M V U = [-M D M^ l M T D ] 



a i 



i=l,2,3 



(6) 



indicating that each sterile neutrino makes a contribution to the active neutrino masses 
of order 8fMi 48|. For example, for an eV-scale sterile neutrino Mj ~ eV together with 
&i ~ 0.3, the contribution to the active neutrino masses is of order 0.1 eV; for a GeV-scale 
sterile neutrino ~ GeV to give a contribution of the same order its corresponding mixing 
angle should be 0, ~ 10~ 5 . As a general rule, the heavier the right-handed neutrino mass, 
the smaller the active-sterile mixing. 

In general the charged lepton mass matrix may not be diagonal: in that case the total 
3x6 lepton mixing matrix connecting the three left-handed lepton doublets L a = (i/ a , a) T 
(a = e, /i, r) to the six neutrino mass eigenstates is 



U 



1 -BB^V V1 V}BV R 



(7) 



where Ve is defined by M e Mj = V £ diag (|m e | 2 , (m^ 2 , \m n 
standard result V^V V for the 3x3 lepton mixing matrix is obtained. 



\ 2 )vl 



Note that for B = the 



C. keV sterile neutrino WDM 

If one of the above-mentioned sterile neutrinos is located at the keV scale and does not 
decay on cosmic time scales, it could be viewed as a WDM candidate. In realistic sterile 
neutrino WDM models, a specific mechanism for the relic production of sterile neutrinos 
is required. For instance, in the Dodelson-Widrow scenario, i.e. production by neutrino 
oscillations, if one assumes that sterile neutrino WDM with mass M s and mixing 8 S makes 
up all the DM in the Universe, its abundance is given by ^, H-Es} 

^^.fJLu^r. ( 8) 



^3 x 10^ 9 / V3 keV^ 

In this work we do not focus on a specific production mechanism of sterile neutrino WDM, 
but will take Eq. ([8]) as a guideline and demand that the WDM neutrino has a mass of a few 
keV and mixing of order 10~ 4 with the active sector. Our main focus lies on the feasibility 
of accommodating sterile neutrinos in flavor models. 

It should be noticed that such a light sterile neutrino (u s ) results in a contribution 9 2 S M S ~ 
10 -5 eV to the active neutrino masses, which is much smaller than the lower bound from 
oscillations of (9(10 -2 ) eV, and hence can be safely ignored when discussing active neutrino 
masses and mixings. Effectively, one can decouple v s in the seesaw formula, leaving only a 
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5x5 mixing matrix together with 2 massive right-handed neutrinos, and a 3 x 5 mixing 
matrix in Eq. (J7|). We present an explicit model in Sect. IHII in order to realize such an 
effective picture. 



D. Neutrinoless double beta decay 

As already mentioned in the introduction, neutrinos with mass below \q\ ~ 100 MeV 
contribute to the Ouf3(3 process via an effective mass defined by (m ee ) = \ Yli=i ^1i m i I > 
where i runs over all the light neutrino mass eigenstates. On the contrary, for right-handed 
neutrinos with masses much larger than |g|, their effect in Ozv/5/5 is strongly suppressed by 
the inverse of their mass. Therefore, if all the right-handed neutrinos are light, i.e. M 2 <C q 2 , 
one obtains 



(m ee ) 



E + E U < 



[M 6X6] 



(9) 



showing that the effective mass cancels exactly, since the the (1, 1) entry of the full 6x6 
neutrino mass matrix in Eq. ([T|) is vanishing. However, this cancellation is not realized if one 
of the right-handed neutrinos is very heavy, since one should decouple this heavy neutrino 
in computing the amplitude for 0z^/3/3. 

The result in Eq. (Q holds in the general framework of type I seesaw models. However, 
in certain flavor symmetric seesaw models in which neutrino mixing is entirely determined 
by the Dirac mass term, Mp can be expressed as 



5J,|55| 



M 



D 



V v diag (y^mM, y/-m 2 M 2 , ^-m 3 M 3 ^ V 1 



T 
R 



(10) 



The active-sterile mixing in Eq. (TJJ is now given by 6 ai = U a ,3+i = (V u ) a i'\/— Ttj/Mj, which 
is merely a rescaling of each column of V„, indicating a direct connection between active and 
sterile sectors. Interestingly, this implies that the above-mentioned cancellation for light 
right-handed neutrinos in (m ee ) occurs pairwise, since 



-(K 2 



1,2,3) 



(11) 



neglecting terms of order B 2 in Eq. ([7]). Here we have assumed Mi to be diagonal, but the 
result still holds with non-trivial Vg, which can be factored out from both U e i and £4,3+4 • Put 
into words, this result means that the contribution to (m ee ) from the z-th active neutrino is 
exactly cancelled by the contribution from the i-th sterile neutrino. This actually simplifies 
the computation of (m ee ) since in Eq. one only needs to count the effects of those active 
neutrinos whose corresponding sterile neutrinos are heavier than \q\. 
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III. At SEESAW MODEL WITH ONE keV STERILE NEUTRINO 



In this section we describe an seesaw model with three right-handed neutrinos: one 
at the keV scale and the other two at either the eV scale, the heavy scale GeV), or both. 
The FN mechanism is used to control the mass spectrum of right-handed neutrinos and 
to set the charged lepton mass hierarchy; since most A4 seesaw models place right-handed 
neutrinos in the triplet representation (see the classification table in Ref. 27|) one has to 
make non-trivial modifications to those models in order to assign different FN charges to 
each sterile neutrino 4 . Indeed, in order to get TBM 57| . 



/ 



U' 



TBM 







\ 



V 



(12) 



2 / 



at leading order with diagonal right-handed neutrinos as A4 singlets, one must choose the 
vacuum expectation value (VEV) alignments of the flavon fields along the directions of 
the columns of the TBM matrix, similar to the method outlined in Refs. 58M60J. The 
crucial point is that each light neutrino mass eigenvalue rrti is then suppressed by only one 
of the heavy right-handed neutrinos M i; so that one can decouple any one of the right- 
handed neutrinos and still achieve TBM with the remaining two columns, at the price of 
one massless active neutrino. Since W12 7^ 0, it is only viable to decouple the neutrinos 
that correspond to the first or third columns, giving normal {m\ = 0) or inverted (m^ = 0) 
ordering, respectively. The decoupled right-handed neutrino becomes the WDM candidate. 

In what follows, we will show a concrete model example in the type I seesaw framework, 
and outline various possible scenarios that differ by the mass spectra of both active and 
sterile neutrinos. In each case we demand one right-handed neutrino to be at the keV scale, 
whereas the other two could be at very different scales, depending on the chosen FN charges. 
Each scheme exhibits distinct phenomenological signatures. 



A. Outline of the leading order model 

Here we outline the model and give general analytical results, focussing on the decoupling 
of the WDM sterile neutrino. Tabled shows the particle assignments of the A 4 seesaw model, 
with right-handed neutrinos uf (i = 1, 2, 3) transforming as singlets under A4. Three triplet 
flavons ip, ip' and cp" are needed to construct the columns of Mp as well as the charged lepton 
mass matrix, and the singlet flavons £, £' and £" are introduced in order to give masses to 



561 ] also has right-handed neutrinos as singlets, but instead of the FN mechanism a 



The model in Ref. 
hierarchy amongst the flavons is assumed. 
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TABLE I: Particle assignments of the A$ type I seesaw model, with three right-handed sterile 
neutrinos. The additional Z 3 symmetry decouples the charged lepton and neutrino sectors; the 
£/(1)fn charge generates the hierarchy of charged lepton masses and regulates the mass scales of 
the sterile states. 



Field 
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e c 


If 




h u ,d 




</ 




e 


v\ 


" C 2 




SU(2) L 
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1 


2 


1 1 
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1 1 1 
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1 
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1 


1" 


1' 
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3 3 


3 


1 i' 1 


1 


1 


V 


1 


z 3 


UJ 


uj 2 


uj 2 


UJ 2 


1 


1 UJ 


UJ 2 


CO 2 UJ 1 


1 


uj 2 


UJ 


1 


tf(l)FN 




3 


1 













-1 


Fx 


F 2 


F 3 



the right-handed neutrinos and keep Mr diagonal at leading order. The NLO terms implied 
by the presence of the flavons will be discussed later. The Lagrangian invariant under the 
SM gauge group and the additional At (g) Z 3 ® £7(1 ) FN symmetry is 



A -A 3 (tpLh d ) e c + ^A (tpLhi)' n c + ^ (tpLh d )" r c 
V j-\ F ^Lh u )v{ + V j-\ F *{<p'Lh u )"v c 2 + V j-\ F *{<p"Lh u ) 



(13) 



at leading order, where the notation (afe)' refers to the product of A 4 triplets transforming 
as V, etc., and y^, and are coupling constants. A = (6}/A < 1 is the FN suppression 
parameter, and for simplicity we assume A to be the cutoff scale of both the v4 4 symmetry 
and the FN mechanism. 

If one chooses the vacuum alignment 5 (<p) = (v,0, 0), the charged lepton mass matrix is 
diagonal 6 : 



fyA 3 \ 



VdV 

A 



V 



^A 
y T J 



(14) 



where v d = (hd) and the charged lepton mass hierarchy is generated by the FN mechanism. 
The right-handed charged leptons e c , /i c and r c carry different charges under the £/(l) FN 
symmetry (cf. Tabled]), which leads to their observed hierarchy. We will employ the same 



5 Note that our model contains two Higgs doublets for the up- and down-sector, respectively, and therefore 
can be accomodated within supcrsymmctry. The VEV alignment could in this case be arranged by 
"driving fields" 61 1. 

6 NLO operators will modify the structure of Mi, introducing non-trivial mixing in the charged lepton 
sector (see the Appendix). 
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mechanism in the right-handed neutrino sector; for the moment the FN charges of the right- 
handed sterile neutrinos are left as free parameters, allowing us to discuss different mass 
spectra. 

As discussed in Sect. [Til a sterile neutrino with mass M, = 0(keV) and mixing of order 
Of ~ 10~ 8 will give a negligible contribution to neutrino mass, and can thus be decoupled 
from the seesaw mechanism. It is then expedient to work in a 5 x 5 basis, with the Dirac 
mass matrix Mo a 3 x 2 matrix and Mr a 2 x 2 symmetric matrix. This is analogous to 



the minimal seesaw model |62l . |63| and the z/MSM, in which the lightest active neutrino is 
massless. The mass spectrum of active neutrinos can either have normal ordering (NO), with 
m 3 ^> m 2 ^> rrii — 0, or inverted ordering (10), with m 2 ^ ra\ ^> m 3 ~ 0. However, there 
exist different scenarios depending on the FN charges assigned to the remaining right-handed 
neutrinos. In order to keep the presentation concise we give general analytical formulae in 
this subsection and discuss details specific to the mass spectrum later on. 

In our model, v\ is assumed to be the WDM candidate, with a mass given by 

Mi = Wl u\ 2F \ (15) 

where u = (£). Note here that Majorana mass terms are doubly suppressed by the FN 
charge. The vacuum alignment (if) = (v, 0, 0) means that at leading order the first column 
of the Dirac mass matrix in Eq. ([Q) is (yivv u X Fl / A, 0, 0) T , so that the sterile neutrino 
only mixes with the electron neutrino 7 . From Eqs. (jlj) and ([5]), the active-sterile mixing is 

0* - = ^A-* , (16) 

so that the FN charge F\ actually enhances the active-sterile mixing, and the contribution 
of the sterile neutrino v\ to the lightest neutrino mass is 

ml < 3 = ^- (17) 

Once we fix the scale of the various flavon VEVs, Fi is fixed by the WDM constraints [which 
we assume to be the ones in Eq. (JSJ)], and the various scenarios to be discussed will differ 
only by the choice of the FN charges F2 and F3, i.e. the scale of the remaining two sterile 
neutrinos. 

With the keV sterile neutrino v\ decoupled, the seesaw proceeds with the remaining two 
right-handed neutrinos, z/£ and v^. For the NO case, we assume the triplet VEV alignments 8 

(<^') = (<;»'), <¥>"> = (o,«", -tO , (is) 



NLO 
Ref. 



terms will induce mixing between v\ and v^ T (cf. Sect. UTTB)) . 



56j employs a radiative symmetry breaking mechanism in order to achieve this VEV alignment. 
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which result in the following 5x5 neutrino mass matrix in the basis (z/ e , v^, u T , v^, 



m5x5=( M D 
\Ml M R 



(19) 



with the Dirac mass matrix 



M 



(NO) 
D 



Vu 

A 



fy 2 v'X F2 \ 

y 2 v'X F2 -y 3 v"X F * 
\y2v'X F2 y 3 v"X Fs ) 



(20) 



and the right-handed neutrino mass matrix 



w 2 u'X 2F2 







w 3 u 



»\2F 3 



(21) 



where u' = (£') and u" = (£"). 

The neutrino masses and flavor mixing can be obtained by the full diagonalization of 
M^ x5 , i.e. UlM^ x5 U* = diag(m 1; m 2 , m 3 , m 4 , m 5 ), where m 4 and m 5 denote the masses of 
right-handed neutrinos. Since eV-scale sterile neutrinos may be present, one should include 



the N 
Refs. 



O seesaw terms, as motivated above. Using the formalism outlined in Eq. (j2J) and 
4451]. and assuming real matrices for simplicity, one arrives up to order ef at 
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-4) 



where the expansion parameters are given by 



Cl = ^%A-^ and 
w 2 u'A 



y 3 v"v. 
£2 = - 



w 3 u"A 



(22) 



(23) 



in analogy to Eq. (j4]). These parameters control the size of active-sterile mixing and NLO 
corrections to neutrino masses and mixing, and will be important in the discussions of various 
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scenarios in the following subsections. The neutrino mass eigenvalues are 



TTh\ 


- n 

— u > 




m 2 


= mf (1 - 


3e?) , 


m 3 


= mf (1 - 


2e 2 2 ) , 


mi 


= w 2 u'\ 2F * 


(0) 

— m 2 


m 5 


= w 3 u'\ 2F * 


(0) 

— m 3 



(1 - 3el) 
(1 - 2ej) 



plus higher-order terms, where 



m 



(0) 



w 2 u'A? 



ni 



(0) 



are the leading order seesaw terms in the NO. 

For the 10 case the following VEV alignments are assumed: 

(</) = <</) = (2<A -t 

The Dirac mass matrix is modified to 



-v") . 



M (io) = ^ 



A 



/ y 2 v'\ F2 2y 3 v"\ F A 

V2v'\ F * -y 3 v"\ F * 
\y 2 v'\ F2 -y 3 v"\ F *J 



(24) 



(25) 



(26) 



(27) 



while the right-handed neutrino mass matrix Mr remains unchanged. In this case, the 
diagonalization matrix approximates (up to order e 2 ) to 
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(28) 
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and the neutrino masses are given by 



mi 


= mi 0) (1 - 

1 \ 


Gel) , 

2 ) ' 


m 2 


= mf (1 - 


36?), 


m 3 


= 0, 




m 4 


= w 2 u'\ 2F * 


(0) 

— m 2 


m 5 


= w 3 u"X 2Fs 


(0) 

— m\ 



(29) 

(1 - 36?) , 
(1 - 6e 2 2 ) , 



(30) 



where 

m (o)_ Gyjv'-vl 

1 — o 
W? l u"h jl 

is the leading order expression for the lightest mass in the 10 and m 2 is defined in Eq. (|25|) . 
Note from Eqs. ( 122]) and ( 128]) that the mixing pattern \U e %\ = and (L^l = \U T s\ is stable 
with respect to higher order seesaw terms, which is actually true to all orders in e, 45 1. 

One salient feature of the above seesaw model can be seen from Eqs. ( 125]) and ( 130]) : the 
leading order contributions to the active neutrino masses do not depend on the FN charges 
assigned to the right-handed neutrinos. The leading seesaw mass term is M^/Mj, so that the 
one unit of FN charge X Fi from Md cancels with the two units X 2Fi from Mj. On the other 
hand, the NLO term M^/Mf oc e 2 does depend on the FN charge, which therefore controls 
the magnitudes of NLO corrections. The larger the charge (equivalent to a smaller sterile 
neutrino mass), the larger the correction parameters q become, and thus the larger the 
corrections to the leading order seesaw masses. 

In addition to NLO seesaw terms, one would expect higher-dimensional operators to 
modify the leading order predictions of the model, which has so far been constructed from 
the leading order Lagrangian in Eq. (TT3"]) . The magnitude of those corrections depends 
largely on the actual numerical values chosen in the model, since they are suppressed by 
additional powers of the cutoff scale A. Our choice of mass scales is guided by the leading 
order predictions: we need (i) the sterile neutrino mass and mixing to satisfy Eq. ([8]), (ii) the 
correct scale of active neutrino masses and {Hi) Yukawa couplings to be < 0(1). In what 
regards the keV sterile neutrino [see Eqs. ( fT5l) and ( fl6l) ]. a rough numerical estimate shows 
that with the mass scales 

v ~ 10 11 GeV, u ~ 10 12 GeV, A ~ 10 13 GeV , (31) 

the Higgs VEV v u = (h u ) ~ 174 GeV and A ~ 0.1, one needs the FN charge 

Fi = 9 (32) 

to obtain a sterile neutrino of mass Mi ~ 1 keV with the desired mixing angle Q\ ~ 10 -8 , 
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with yi, w\ < 0(1). In order to stabilize the active neutrino masses around the sub-eV scale, 
one can choose [together with the numbers in Eq. fl3"Tj) ] the scales 

v' ~ v" ~ v! ~ u" ~ 10 11 GeV (33) 

for the flavon VEVs, and the mass splitting among active neutrinos can be achieved by prop- 
erly choosing the corresponding Yukawa couplings, i.e. and Wi (i = 2, 3). For definiteness 
we fix the scales of the VEVs from here on, and obtain all numerical results using those 
values. 



B. Mixing corrections from higher-order terms 

As we have already mentioned, the presence of gauge singlet flavons in the model will 
inevitably induce NLO corrections, which may modify the leading order picture and affect 
both active and active-sterile neutrino mixing. Indeed, modifications to TBM are required 
in order to explain the T2K result that suggests non-zero #13 64J. We concentrate on 



the effects of adding higher-order operators to the Lagrangian in Eq. 
introduce corrections by perturbing the A4 triplet VEV alignments 
getting non-zero #13 in models designed to predict TBM is a more general problem, and 



TT3"]) ; one could also 
65j. Note that 



27 



other solutions have been proposed, e.g. in Refs. [66 



681. 



Since the charged lepton and right-handed neutrino mass matrices are diagonal at leading 
order, TBM comes solely from the structure of the Dirac mass matrix. Without performing 
a detailed numerical analysis, one can show that the higher-order corrections affect all three 
mass matrices: Me, Mb and Mr. The impact of those corrections is controlled by the ratios 
of flavon VEVs to the cut-off scale, in our case 



U U U V V 

r\ = — ~ 0.1 and r 2 = — — — — — — — — 0.01 . 



A 



A A A A 



(34) 



The terms containing the VEV (£) = u = r%A have the largest effect, and will be included 
in our analysis (see the Appendix); terms containing the VEVs it', u" , v, v' and v" are all 
of relative order r 2 — 0.01 and can be safely neglected. Importantly for our model, the 
correction terms turn out to have a negligible effect on the keV sterile neutrino mass, as well 
as its mixing with the active sector. Explicitly, from Eqs. (IA-121) and (IA-141) . the corrected 
active-sterile mixing is 



0' 



(NO) 



el 



6el 1 



y_W 
yiv 



-n 



and Q' el 



(IO) 



del 



yW 



2 y J y< ] 

yiv wi 



(35) 



where the dimensionless couplings y[ and w[ are defined in Eqs. flA-51) and flA-81) . respectively, 
and the leading order expression for # el is given in Eq. ( I16p . In addition, the mixing angles 
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Bnx and 6 Tl become non-zero, but of the same magnitude as 9 e i, i.e. 



a> (NO) 



del 



yiv 



yiv w 1 



n and 6[ lTl m 



yW 
yiv 



yiv wi 



(36) 



This shows that the active-sterile mixing is stable, illustrating the point that unlike active 
neutrino mixing it is defined as the ratio of two large scales, so that small changes in Mu 
and Mr will have little effect on 9 a i (we assume that \w' x \ ^ | wi | ) . The WDM particle 
remains decoupled from the seesaw and one can still work in the 5x5 basis. We show the 
resulting mixing matrix elements here and provide details of the diagonalization procedure 
and modified neutrino mass eigenvalues in the Appendix. 

The final lepton mixing matrix is a 3 x 5 matrix connecting the three flavors of lepton 
doublets to the five neutrino mass eigenstates, and corrections from the charged lepton sector 
[Eq. f ]A-3j) ] and the neutrino sector [Eq. ( 1A-15|) ] can be combined via Eq. ([7]) to give the 
approximate mixing matrix elements 



\u, 



(:2\ 



\u, 



fj,3 I 



\Ue M \ 2 ^ 4 

\U J 2 ~ f 2 

\ u e5\ — fc 2 



Ij'r 



y^ y-r 



+ l(x-p 3 ) 2 -(x-P3)n C^-* 
2 \y» y T 



2e 



y T 
y T 



2 

3 



(Hk. 


.£)} 




yrj . 


T" 


1 



1 T 2p2 T 2n 



\U, 



fj,5 



2 ^ J2 

Co 



y» 
o y'r 

LT\ — 

yr 



Un y 

2 

2ri 1 



it 



{x - Ps) + (x - Psf 



(37) 



in the NO. Here the are generated by NLO seesaw terms, y'^ T stem from corrections to 
the charged lepton mass matrix, while the other parameters come from corrections to Md 
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and Mr. For the inverted ordering we find 



\u. 



c2 I 



If/, 



f*3| 



It/, 



c5 



1^ 12 



V5| 



2 
1 
3 
1 
2 



4e 



1 - 3e 2 - 2p 2 
1 + 2p 2 
1 - 2p 2 =f 

l + n 
- 2n 



P2^1 — 

Mil Vr 

V» Vr 



pi 



2 ^r 



2 — T\ 

Vr 



(38) 



Vfi Vr 
I 



Vr 

y 



(x - Pa) 



with the parameters 



CT± = X ± P2 ~ P3 , cri = X ± 3p 2 - P3 



X 



yiv w 1 

y 3 v" wi 



-n , P2 



R = 



in 



(o) 



(0) 



Am| 



2/2^ 

T r i' 

c^io -1 ; 



P3 



(39) 



G = 



m 



(o) 



m 



(0) 



(o) 



2Ami 
Arrtg 



i? 2 



C(10 5 



controlling the size of the mixing terms, where Am 2 , and Am^ are the solar and atmo- 
spheric mass squared differences, respectively. The dimensionless couplings y' 2 3 are defined 
in Eq. (1A-5|) . R and G contain the leading order neutrino masses from Eqs. (125]) and (130]) : 
while R is quite small, G is large, which is a consequence of the two relatively large but 



nearly degenerate neutrino masses in the IO, m\ 



(0) 



m 



(0) 



0.05 eV. We have expanded to 



first order in R, but G remains an exact expression in the mixing matrix. Thus in the IO 
we need cr^ = x + 3p 2 _ P3 to be C(10^ 3 ) in order to keep the corrections to \U e 2\ 2 under 
control, which in turn puts a constraint on the Yukawa couplings y' 2 3 and w[ in Eqs. ( 1A-5I) 
and (IA-81) . With n ~ 0.1 and u ~ ?/ ~ v", we have P23 — 0.1— and y ~ 0.1^ — , so that 

' 1 ' 2/2,3 A !/3 Wi ' 

we need to assume that y 23 — 0.01?/ 2j 3 and yin^ — O.OlysWi in the inverted ordering. The 
full neutrino mass eigenvalues are given in Eqs. ( ]A-17j) and (1 A- 19ft : despite the appearance 
of G 2 terms in the IO mass eigenvalues they will always be suppressed by (crj_) 2 , which is 
constrained to be small from the mixing matrix element U e 2- 

As expected, by setting y' 2 , y' 3 , w[, y'^ and y' T to zero in Eqs. ( 1371) and ( 1381) one recovers 
the matrix elements in Eqs. (122]) and (128]) . Note that without the higher-order correction 
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terms U e3 remains exactly zero, to all orders in The active-sterile mixing (U a ±$) is always 
proportional to 6j, or in other words to a ratio of scales [cf. Eqs. (jl]) and (|23j) ] . In the different 
scenarios discussed in the following subsections, the terms will have different magnitudes, 
depending on the right-handed neutrino spectrum. In those cases with significant values of 
€i (i.e. eV-scale sterile neutrinos) one must take into account both NLO seesaw corrections 
and higher-order corrections, whereas in cases with negligible e, (heavy sterile neutrinos) 
one need only worry about the higher-order correction terms, i.e. those controlled by y' 2 , y 3 , 
w[, y'^ and y' T . 

Even if y' 2 , y' 3 and w[ are small and mixing corrections from the neutrino sector are 
negligible, there are still effects from the charged lepton sector. Indeed, in order to keep 
the solar mixing angle within its allowed range 0|, one has the constraint (assuming for 
definiteness y' 2 = y 3 = w[ = and — 0) 

_0.4< ( V -» + y± \ < 0.95, (40) 

on the charged lepton Yukawa couplings; the extreme choice y'Jy^ = —y' T /y T gives the 
reactor mixing angle 

sin 2 # 13 ~ 1r\ (^\ ~ 0.02, (41) 

in both mass orderings, assuming that y' T ~ y T . In this case sin 2 # 2 3 — 0.6, and sin 2 # 12 retains 
its TBM value. 



C. Explicit seesaw model scenarios 

In order to illustrate the versatility of the model discussed, we present three scenarios 
with different mass spectra in the right-handed neutrino sector. Each case differs by the 
choice of FN charges F2 and F 3 , what one could call the "theoretical input"; the consequent 
neutrino phenomenology is described in detail. Table [Til summarizes the key differences in 
each case. 

In all cases we have checked that Yukawa couplings of order 1 or 0.1 can fit the model 



to the active neutrino mass-squared differences 47]], and, where appropriate, to sterile mass 



parameters 30|. The effects of the higher-order corrections discussed in Sect. IIII Bl are 
described for each scenario. Due to the large number of parameters we will always have 
enough freedom to fit the masses to the data, so that we only need to take care that mixing 
corrections are under control, particularly in the IO, as discussed above. 
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TABLE II: Summary of the different scenarios discussed in the A4 seesaw model. In each case 
the WDM sterile neutrino has a mass Mi = O(keV), and the corresponding active neutrino is 
approximately massless. 





\7Tl ) 

F t , F 2 , F 3 Mass spectrum \U a A \U a5 \ x ee/ 

NO 10 


Phenomenology 


T 
1 


9, ID, ID iV22,3 = C(ev), C(O.l) C(O.l) U 


3 + 2 mixing 


IIA 


9 in n M 2 = 0(eV) romiWonn -n^ n 2y/Am\ 
M 3 = O(10 n GeV) 3 


3 + 1 mixing 


IIB 


qnm M 3 = 0(10" GeV) „ rin _^ n(n „ ^| V / A ^i 
M 3 = O(eV) 3 3 


III 


9,5,5 M 23 = O(10 GeV) C(10~ 6 ) C(10~ 6 ) ^Ami 


Leptogenesis 



1. Scenario I: two eV-scale right-handed neutrinos 

In this case we assign the FN charges F± — 9, ^3 = 10, so that the right-handed 
neutrino masses are lowered down to the eV scale. It is now notable that e\ )2 = O(0. 1) 
can be expected, indicating that NLO seesaw terms should be considered. The effects are 
more pronounced in the IO case, since two of the active neutrinos are nearly degenerate and 
are more sensitive to corrections. The five neutrino mass eigenvalues are given by the full 
expressions in Eqs. ( 1A-17I) and ( 1A-19I) . 



In this scenario, there are no heavy right-handed neutrinos that could be used to ex- 
plain the matter-antimatter asymmetry via leptogenesis. Neutrinoless double beta decay is 
also vanishing since the contributions from active and sterile neutrinos exactly cancel each 
other 9 , unless there are other new physics contributions. However, the eV-scale right-handed 
neutrinos offer an explanation for the short-baseline oscillation anomalies often attributed 
to them. 

In the NO case, one of the two sterile neutrinos could mix with v e via ?7 e4 ~ e±. The 
reactor flux loss is therefore explained since part of the total flux of V e oscillates into sterile 
neutrinos. However, one finds that the active-sterile mixing turns out to be too tiny to 
account for the reactor anomaly. This can be deduced from Eqs. (I24p and (129]) : at leading 
order, e\ oc m 2 /m 4 . In the NO, m 2 — 0.009 eV is fixed by the neutrino mass-squared 
differences, and hence, e% can hardly be sizable for an eV-scale m 4 . The situation is different 
for the IO case, since m 4 ~ m 2 ~ 0.05 eV is fixed from neutrino oscillation experiments. 



9 This is different to the usual analysis, e.g. in Ref. 25( (see also 69|), in which sterile singlet states are 
simply added to an existing model. 
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FIG. 1: The allowed ranges of \U e ^\ 2 — Am^ (blue) and \U e ^\ 2 — Am^ (red) in the inverted ordering, 
requiring that the oscillation parameters lie in their currently allowed 2a ranges. The blue and 
red vertical and horizontal error bars indicate the allowed 2a range for the 3 + 2 mass and mixing 
parameters from Ref. 13(1 . their intersection is the best-fit point. The black errors bars are for the 
3 + 1 case from Ref. [30|, to be discussed in scenario II in Sect. IIIIC 21 



Furthermore, both U e 4 and U e ^ are non- vanishing (see Fig. [I]). 

The effect of higher-order operators on the active-sterile mixing is very small. Switching 
on w[ gives \U e ^\ 2 ~ 0(r\) e\ in the NO [cf. Eq. fl37|) ]. which will still not give sufficient 
mixing to explain the data. In the 10 case, |t/ e5 | 2 ~ 4[1 + 0{ri)]e\, so the small correction 
term makes little difference. Indeed, the allowed ranges illustrated in Fig. [1] already include 
the effects of higher-order operators. One observes that the desired active-sterile mixing can 
indeed be achieved in the 10 case. 

In what regards active neutrino mixing, deviations from TBM come from both NLO 
seesaw terms (oc 6j) and higher-order operators (oc y 2 ,y' 3 ,w' 1 ,y' ll ,y' T ). If we only consider 
higher-order corrections in the neutrino sector for simplicity, i.e. the y' 2 3 and w[ terms in 
Eqs. (IA-5j) and flA-8j) respectively, then from Eqs. f )37|) and f )38|) only oc e 2 receives 
visible corrections in the NO, since e\ and the product a+R are both small. However, the 
higher-order terms related to the product a+G lead to sizable corrections to \U e 2\ 2 in the 
10 case; \U e2 \ could be enhanced or reduced depending on the signs and magnitudes of y' 2 3 
and w[. In addition, non-zero #13 can be obtained from the charged lepton corrections, as 
discussed in Sect. IIIIBI above. 



20 



2. Scenario II: split seesaw with both eV-scale and heavy right-handed neutrinos 



We have shown that it is possible to get either normal or inverted ordering by choosing 
the alignment of the flavon VEV (<//') correctly [cf. Eqs. (jl8p and (12 6p ]. In this case we now 
assign different FN charges to the two seesaw right-handed neutrinos, so that there are four 
distinct possibilities, depending on the mass ordering of active neutrinos and which sterile 
neutrino {v 2 or h>Q is chosen as the heavy one. One can then use a two-stage seesaw, by 
integrating out the heavy sterile neutrino first and then applying the seesaw formula again. 
With the assignments F\ = 9, F 2 = 10 (0) and F3 = (10) the sterile neutrino z/3 (z/|) 
has a mass in the 10 11 GeV range, and is integrated out first, whereas v 2 (^3) is at the eV 
scale. The third (second) column of Mr> is then used in the seesaw formula, leading to a 
3x3 effective neutrino mass matrix of rank 1 that gives one of the active neutrinos masses. 
The full 4x4 mass matrix in the basis (u e , u T , ff/ 3 leads to mixing between the active 
sector and the remaining eV-scale sterile neutrino v 2 (^3)- Here one can apply the method 
and formulae outlined in Sect. IIII Al except that one has a 4 x 4 mixing matrix, which can 
simply be obtained from the formulae in Eqs. ( 1221) and ( 1281) by removing the relevant row 
and column. 

• Case IIA: v\ heavy, (F 3 = 0), v c 2 light (F 2 = 10) 

In this case one removes the fifth row and fifth column of U v in Eqs. ( 1221) and (1281) . 
giving the same 4x4 mixing matrix in both mass orderings, and the matrix elements 
\U e §\ 2 and It^l 2 are zero. The light neutrino mass eigenvalues (i = 1,2,3,4) are 
given by the expressions in Eqs. (1A- 17[) and (1A- 19[) with e 2 set to zero; the heavy 



neutrino has the mass M3 = w^u" . It is the small value of F3 that leads to e 2 — 
[Eq. ( |23l) ]. so that (or mi) does not receive any higher-order corrections, as this 
mass originates from the high-scale of M 3 , whose FN charge "cancelled" in the leading 
order seesaw formula. Although in our example we have F3 = 0, so that M3 ~ 10 11 
GeV, even with F3 = 5 and M3 ~ 10 GeV, one has e 2 — 10~ 6 (see scenario III), so 
that NLO corrections would still be under control. 

The FN charge F 2 = 10 of the eV-scale neutrino gives corrections to m 2 and m^, via 
ei ~ 0.1. With order one Yukawas and values for the VEVs as before, M3 lies around 
10 11 GeV. The effective mass in 0^/3/3 is given by the (1, 1) element of the 4x4 mass 
matrix, which, at leading order, is 



2mS 0) 



2 

m 



0.032 eV. (42) 



(m ee )( N °)=0, (m ee )^ 



Here one can see that the contribution of the light neutrino of mass m 2 has cancelled 
with that of the light sterile neutrino z/g, in both mass orderings. Note again that 
this is different from the usually discussed effects of sterile neutrinos in The 
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effective mass is zero in the NO since at leading order, U e s = 0. A non-zero value 
of U e s would give a very small contribution to the effective mass in the NO, and a 
completely negligible one in the 10. 

Case IIB: v% heavy (F 2 = 0), i/| light (F 3 = 10) 

Here the mixing matrix is found by removing the fourth row and column of Eqs. (F22J) 
and ( 128]) . so that the matrix elements in Eqs. (|37|) and ( 138]) can be relabelled \U e ^,\ 2 — > 
| t/e4 1 2 and |t/ M5 | 2 — > \U^\ 2 . The light neutrino mass eigenvalues m; {i = 1,2,3,4) 
are now found by setting e\ to zero in Eqs. (IA-17P and (IA-191) . with the relabelling 
m 5 — y m 4 ; the heavy neutrino has the mass M 2 = W2U'. The roles of the sterile 
neutrinos are now swapped, and M 2 is situated at the 10 11 GeV scale. The effective 
mass at leading order is 



(m ee > (NO) 
<m ee > (IO) 



in 



(0) 



m 



(0) 



0.0029 eV 



0.016 eV 



(43) 



in this case the contribution of has cancelled. Again, corrections to the mixing 
angles give very small corrections to the effective mass. 

In both cases IIA and IIB one could potentially explain the reactor anomaly in the 
framework of 3 + 1 neutrino mixing 30|, |70[, with \U e 4\ 2 ~ [1 + O(^i)] e\ in case IIA and 



\U e 4\ ~ 4[1 + 0(ri)]e 2 in the IO in case B. Once again only the IO fits the data: the allowed 
ranges in the mass-mixing plane for the IO in case IIA (IIB) are shown by the blue (red) 
points in Fig. [TJ One can see that the best-fit point (the black cross) from Ref. |30( is 
compatible with case IIB. Finally, the effects of higher-order operators on both active-sterile 
mixing and active mixing are the same as in scenario I, except that one should switch off 
the effect of e 2 (ei) in case IIA (IIB). 



3. Scenario III: two heavy right-handed neutrinos 

In this case we take iq = 9, F 2 ^ = 5, so that one can estimate that the q ~ 10~ 6 
(i = 1, 2) are dramatically suppressed, and the NLO seesaw terms in Eqs. (122]) and (128]) can 
be safely neglected. The 3x3 effective neutrino mass matrix is given by Eq. (J3j), with 
defined in Eqs. ( [201 or (1271) and Mr from Eq. ([21]) ; the active neutrino masses are simply 
given by the leading order masses mf\ The heavy neutrinos have masses M 2 = w 2 u'X 10 
and M3 = w 3 u"\ w . Without the effect of the terms, the only modifications to the TBM 
pattern come from the higher-order operators in Sect. IIII Bl 
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The two heavy right-handed neutrinos that participate in the seesaw formula have masses 
around 5 GeV, assuming order one Yukawas and the usual values of the VEVs. Note that 
one could set W2 = w% to obtain degenerate right-handed neutrinos M 2 = M 3 . The choice of 
degenerate sterile neutrinos in the few GeV regime would correspond to the i^MSM paradigm, 
in which no new scales between the SM and the Planck scale are assumed. Baryogenesis 
then proceeds via oscillations between and 1/3, which need to be sufficiently degenerate 
\\M-2 — M 3 |/M 2 ~ 10~ 6 ) to give the correct baryon asymmetry [7l| . 

If we choose = instead, then M.2$ — 10 11 GeV, so that the CP-violating decay 
of right-handed neutrinos could explain the matter dominated Universe via thermal lep- 
togenesis. The required CP violation may originate from complex Yukawa couplings. We 
further note that, similar to the ordinary type I seesaw, neutrinoless double beta decay is 
allowed, and the right-handed neutrinos play no role in this process since their contribution 
^2 i=2 3 &i/Mi is strongly suppressed by the inverse of their mass. Explicitly, at leading order 
the effective mass from the (1,1) entry of Eq. is 



(m ee )( N °) 



m 



(0) 



mi 



2ml 



(0) 



+ 



(o) 
m\ 



0.0029 eV 



Am 2 , 



0.049 eV 



(44) 
(45) 



where the mass eigenvalues are real. If mf and are complex, the IO case becomes 
(m ee )( 10 ) < \J Am\. Corrections from higher order terms are again small. 



(0) 



IV. AN EFFECTIVE THEORY APPROACH 



In this section we recast the idea presented in Ref. [25[, this time in the context of keV 
sterile neutrino WDM rather than eV-scale sterile neutrinos. A popular flavor symmetry 
model, which predicts TBM and is based on the group A4, is modified in order to accom- 
modate a keV sterile neutrino. Unlike the seesaw model, neutrinos get mass from effective 
operators and only one sterile state is introduced. We also extend the discussion to include 
the effects of higher-order operators. 



A. A4 symmetry with one keV sterile neutrino 



The Altarelli-Feruglio (AF) A 4 neutrino mass model [72J is well known, and at leading 
order gives exact TBM for the lepton flavor mixing matrix. The original AF model includes 
three sets of flavon fields <p, tp' and £ in addition to the SM particle content. We add an 
additional sterile neutrino transforming as a singlet under A 4 and Z 3 , with the £7(1 )fn charge 
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TABLE III: Particle assignments of the A4 model, modified from Ref. 72 1 to include a sterile 
neutrino v s . The additional Z3 symmetry decouples the charged lepton and neutrino sectors; the 
FN charge generates the hierarchy of charged lepton masses and regulates the scale of the 
sterile state. 
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of F s = 8. The relevant particle assignments are summarized in Table IHIl 

As discussed in Ref. [25( , at leading order the Yukawa couplings for the lepton sector read 



-C Y = |A 4 (ipLh d ) e c + |A 2 (fpLhi)' fS + | (<pLh d )"r c + ^(Lh u Lh u ) + ^{ip'Lh u Lh u ) 
+ ^\^{p'Lh u )u s + ^\\ip'ip'Lh u )v s + m a X 16 v°X + h.c, (46) 

where m s is a bare Majorana mass. Note that the invariant dimension-5 operator 
-^X 8 (cp' 'Lh u )v s is not invariant under the Z3 symmetry. With the following vacuum align- 
ments (as in the AF model) 

(<p) = (v, 0, 0) , {iff) = (v\ v\ v') , (0 = u , (h u4 ) = v u4 , (47) 

the charged lepton mass matrix is diagonal [cf. Eq. (I14j) ]. and the full 4x4 neutrino mass 
matrix is 

/a + ¥ -I ~i e\ 



M, 



4x4 



3 3 3 

¥ "-I e 

2d 

3 C 



V 



(48) 



where a = 2x a ^r, d = 2x d ^ L and e = v^eA 8 ^^ have dimensions of mass. The first 
three elements of the fourth row of M 4 * 4 are identical because of the VEV alignment {iff) = 
(v',v',v'), which was necessary to generate TBM in the three-neutrino case; this alignment 
combined with the A4 multiplication rules causes the term proportional to Xf in Eq. (146jl to 
vanish. 

If one assumes that a < m s and expands to second order in the small ratio e/m s , the 
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mixing matrix diagonalizing M 4 * 4 in Eq. (1481) is 25] 



U 



with the eigenvalues 



/ -i- 

1 1 1_ 

__L _L 1 

y 







/0 

^ 

^ 

\ o _yje go/ 



3e 2 



3e 2 

mi = a + d , m 2 = a , m 3 = —a + d , m 4 = m s H 

m.. 



(49) 



(50) 



As we will see, the chosen FN charge forces m s to be at the desired keV scale and sets 
the magnitude of active-sterile mixing, e/m s = O(10 -4 ). This means that the "seesaw 
contribution" (oc e 2 /m s ) to m 2 in Eq. (|50|) is negligible. 



B. Estimation of the mass scales and active-sterile mixing 

In order to examine the viability of the model we provide a rough numerical example. As 
discussed in the original AF model {72)], we assume that (z) the Yukawa couplings y and x 
remain in a perturbative regime; (ii) the flavon VEVs are smaller than the cut-off scale and 
(Hi) all flavon VEVs fall in approximately the same range, and obtain the following relation 
constraining the flavon VEVs: 

0.004 < ^ w \ w ^ < 1, (51) 
AAA 

with the cut-off scale A ranging between 10 12 and 10 15 GeV. We would like to suppress the 
mass of the keV neutrino, while at the same time keep its mixing small enough and satisfy 
the conditions in Eq. floTT) . By choosing the FN charge of v s (i.e. F s = 8) and the mass scales 

u ~ v ~ 10 10 GeV , v ~ 10 11 GeV , A ~ 10 12 GeV , 

v u4 ~ 10 2 GeV , (0) ~ 10 11 GeV , (52) 

which means that A = (0)/A ~ 0.1, one obtains 

/ u \ f v u \ 2 ( 10 12 GeV \" 
V10 10 GeV7 V10 2 GeVy \ A / 



A \ 8 / u \/ u' \/ v u \/10 12 GeV\ 2 
kF 1 ; U0 10 GevJ VlO 10 GeV/ UO 2 GevJ V A ) 

with the assumption that the Yukawa couplings x a ,d,e are of order 1. 



(53) 
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The Majorana mass term m s v c s v c s is doubly suppressed by the {7(1)fn charge. There are 
additional terms that can give a contribution to this mass in addition to the bare term. 
From the particle assignments in Table IHH the leading order contribution to m s reads 

(^) A^if ^ A" , (54) 

so that these terms are suppressed by A 16 , and the resulting Majorana mass can be of order 
keV: 

.ID" J (l^Gev) Y— J keV - ^ 

The active-sterile mixing is given by 

e s = — ^l(T 4 , (56) 
m s 

corresponding to sin 2 ^ ~ 10 -8 , in accordance with the astrophysical constraints discussed 
in Sect. [TTJ It should also be noticed that in this model the charged lepton masses are 



v F 

m a = y a v d -X a ~ 10 



so that we get the correct mass spectrum with the FN charges (F a ) of 4, 2 and for e c , // 
and r c , respectively [assuming y a ^ C(l)]- 



C. Higher-order corrections and non-zero 



'13 



One may also wonder if higher-order terms could lead to significant corrections to the 
lepton flavor mixing and neutrino masses so as to generate a non-zero #13, as suggested 
by the T2K experiment. In general, both the neutrino and charged lepton mass matrices 
receive higher-order corrections, suppressed by additional powers of the cutoff scale A; those 
are the only type of corrections that we consider here. 

In the charged lepton sector, the NLO corrections to Mi come from terms like 

^ [y' e X 4 {ipipLha) e c + ^A 2 (^Lh d )' // + y' T {^Lh d )" r c ] , (58) 

which however replicate the leading order patterns, as in the seesaw model (see Ap- 
pendix Ell)- The NLO corrections to Mi can thus be simply absorbed into the coefficients 



y^ 
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As for the sterile neutrino, the NLO corrections to m s are given by 



A 2? ^ A 2 



u 3 3v' 2 u 



16 



(59) 



in this case the contributions in Eq. ( 1591 are of order 10~ 4 keV, and do not affect the scale 
of m s significantly. Note that the term ^A^fi^VWX i s i n principle also allowed, but 
vanishes after A 4 symmetry breaking, just like the Xf term in Eq. (I46p . NLO corrections to 
the e parameter come from terms like 



x 



(60) 



which lead to 



0.01 



10" 



uv 



10 10 GeV) 2 J V10 11 GeV/ V10 2 GeV 



10 12 GeV\ 3 



A 



eV, 



(61) 



indicating again that the active-sterile mixing is hardly affected. 

The higher-order operators contributing to light neutrino masses are of order 1/A 3 . There 
exist only three such terms that cannot be absorbed by a redefinition of the parameters a 
and b [61[, i.e. 



^(ipip')'(Lh u Lh u )" , -^(ipip')" (Lh u Lh u )' , and ^(ipLh u Lh u ) , (62) 



so that the light neutrino mass matrix is modified to 



/a+f 



_d 

3 

2d _ d 
3 " 3 



V 



2d 
3 



• Vi -3^3 



(63) 



where rji = 2x\^^, t] 2 = 2x2^r- an d f]3 = 2x 3 ^^. For 0(1) Yukawa couplings, one can 
estimate that 



T)i ~ 0.01 



10 11 GeV/ V 10 10 GeV J V10 2 GeV 



2/10 12 GeV\ 3 Tr 
I eV 



A 



(64) 



As a result, the NLO terms may lead to visible modifications to the TBM pattern, in 
particular to 6*13, but on the other hand do not entirely spoil the leading order picture, since 
one always has enough parameters to fit the data. Keeping only first order terms in rji, one 
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obtains 



L 1 / \ 1 

m 1 ~ a + - - (r/i + 772) + -773 , 

m 2 ~ a + r/i + 772 , (65) 
m 3 ~ -a + 6 + - (771 + 772) + ^773 , 



together with the mixing angles 

sin 2 6» 13 



8a 2 




K 1 - 


2773 \ 




36 ,/ ' 


K 1 - 


771 - 772 N 




4a y 


-771 = 


0.1a and 



sin 2 £ 12 ~ ± ( 1 - f£ ) , (66) 
sin 2 6»23 



As one numerical example, we take 772 = —771 = 0.1a and 773 = 0.16, and obtain sin" $13 ~ 
0.005, which is compatible with the current global-fit data at 2a CL. In addition, sin 2 6<??, ^ 
0.53 and sin 2 9 12 — 0.31 are predicted, in good agreement with their best-fit values 471. l73|. 



V. CONCLUSION 

The addition of sterile right-handed neutrinos to the SM is a natural way to explain 
light active neutrino masses via the seesaw mechanism. This works even if the scale of 
the sterile neutrinos is not equal to its "natural value" of 10 10 to 10 15 GeV, as long as the 
Dirac mass matrix can also be suppressed such that M^/Mr is small. At the same time, 
several observations point to sterile neutrinos at the keV and eV scales. Therefore we have 
attempted, as a proof of principle, to construct a seesaw model for neutrino mass and lepton 
mixing that can provide a common framework for all these issues. 

Starting from a flavor symmetry model based on the tetrahedral group A4, we described 
different ways to introduce sterile neutrinos, using the seesaw mechanism (and also an ef- 
fective theory approach). In both cases the Froggatt-Nielsen (FN) mechanism is used to 
suppress the masses of the right-handed neutrinos. We stress that its presence in flavor 
symmetry models can be considered necessary in order to generate the observed strong hier- 
archy in the charged lepton sector. In fact, we utilize the very same FN for both the charged 
lepton masses and the right-handed neutrinos. 

In the seesaw model we studied different possible spectra in the sterile sector: once the 
keV WDM neutrino is decoupled one can have the remaining two neutrinos at the eV scale 
or at a high scale (in our example at either 10 GeV or close to the flavor symmetry breaking 
scale of ~ 10 11 GeV). In each case there are distinct phenomenological consequences, both 
for neutrino mass and neutrinoless double beta decay. In particular, NLO corrections to the 
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seesaw formula need to be taken into account when the sterile neutrinos are at the eV scale. 

Motivated by the recent indications for nonzero 6 13 in the T2K experiment, we examined 
the effect of higher-order terms in both the seesaw model and the effective theory. In general 
active neutrino mixing angles will receive corrections of the same order. We highlighted the 
fact that active-sterile mixing is stable in any seesaw model, being defined as the ratio of 
two large scales. 

Although one can explain both eV-scale and keV-scale sterile neutrinos in a single frame- 
work, it is not possible to have viable WDM, eV-scale neutrinos and heavy neutrinos for 
leptogenesis in a model containing three right-handed neutrinos. However, we emphasize the 
point that if one departs from the common theoretical prejudice of right-handed neutrinos 
residing at around the Grand Unification scale, various interesting model building options 
can arise. Further experimental data in the years to come will put the presence of sterile 
neutrinos at the eV and/or keV scale/s to the test, thus determining whether it is indeed a 
useful enterprise to further pursue this avenue of research. 
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Appendix A: Corrections from higher-order operators in the seesaw model 

Here we give details of the procedure followed to calculate corrections to the lepton mixing 
matrix in the presence of higher-order operators, which affect Mn, Mo and Mr. We only 
take into account corrections of relative order r\ c± 0.1 [cf. Eq. ( )34l) ]. Explicit expressions 
for the corrected neutrino mass eigenvalues are also reported. 



1. Charged lepton sector 

The corrections to Mi from dimension-six operators come from coupling a second A4 
triplet or an singlet to each mass term. The addition of the flavon ip replicates the 
leading order pattern, since the triplet from the product (<p<p)3 has a VEV in the same 



61| . Terms with the additional singlet £" also leave the structure of the mass 



direction as <p 
matrix unchanged, but the additional terms 



|fAW^K ,||A 3 £VWe c and V ^\\^ V "Lh d )e c (A-l) 
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are also present, for all three flavors. The first term gives the largest NLO contribution, i.e. 



SM, 



(i) 



VdUV 



y'e^ yft y'r 

\y' e X z y'^X y' T J 



(A-2) 



of relative order n — 0.1. The matrix diagonalizing (Mi + SM^)(Mi + SM^Y can be 
approximated by 
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1 








nk ri 


l 




\- 


fn 

Vt 


1 / 



0(rl\ 2 ) 



(A-3) 



and the charged lepton masses become 



< = (y a + y>i)^V , 



[a = e, /x, rj 



(A-4) 



which amounts to a rescaling of Yukawa couplings. 



2. Neutrino sector 

Similarly to Mg, corrections to Mp from adding the singlet £" retain the leading order 
form, but there are also several terms with two triplet flavons. The latter are all suppressed 
by T2 — 0.01 and can be safely neglected. Of the nine different invariant dimension-six 
operators with one triplet and one singlet flavon, there are three of relative order T\ ~ 0.1, 
namely 



'3 ) 



leading to the corrections 



(A-5) 



SM 



(IN) 
D 



A 2 



y[v' -y' 2 v" y' 3 v 
y[v' y' 2 v" 

\yW 







F and SM 



D 



A 2 



fy[v' -y' 2 v" y' 3 v\ 

y[v' -y' 2 v" 
\y[v' 2y' 2 v" j 



(A-6) 



in the normal and inverted ordering, respectively. Here the matrix of FN charges is 



F = diag(A Fl , X F2 , X Fi ) 



(A-7) 



The corrections to Mr come from terms with two singlets and those with two triplets, 



e.g. 



^ X Fl+F3 ^M + --- and ^A Fl+F3 ( W ')^ 3 c + --- 



(A-8) 
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the singlet terms give the contribution 



/WA 2Fl 



5M R 1] oc - 
A 



V 



u'u"\ 2F * u'u' 



u 



uu\ Fl+Fi \ 

^F 2 +F 3 

"u"\ 2F * J 



(A-9) 



whereas the triplet terms are all suppressed by r 2 ~ 0.01. Comparison of the LO and NLO 
terms shows that the large ratio T\ — 0.1 only occurs in the (1, 3) element of 5M R \ whereas 
the diagonal and (2,3) elements receive small corrections of order r 2 ~ 0.01. Ignoring the 
latter, the new mass matrix is 



M' R = M R + 8M R 1] = F 



1 ' W\u w\ur\\ 



w 2 u 



\ ■ 



(A-10) 



w 3 u" J 



It is convenient to factor out the FN charges here, since they do not appear in the leading 
order seesaw formula. However, as emphasized before, they will play a role when considering 
NLO seesaw terms. Expanding in the small ratios r\ ~ — 0-1, the matrix diagonalizing 
M' R can be approximated as 



V, 



R 



( i o -^n\ 

1 
V^i 1 J 



WiU 



i ) ' 



(A-ll) 



with the mass eigenvalues 



M[ = Wl u\ 2F ^ ( 1 + 
M' = w 2 u'\ 2F2 , 



(A-12) 



ML = W3U 



This shows that corrections to the masses M\^ are suppressed by rf, and the WDM candidate 
v\ remains in the keV range. 

The diagonalization matrix in Eq. flA-llj) can be absorbed into Mp, so that the leading 
order neutrino mass matrix is 



M' u = -M' D &zg{M'-\M'-\M'f x )M' D T 



(A-13) 



where M' D = ( Mjj + 5M^ ) V R and the FN charges have cancelled. The Dirac mass matrices 
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in Eqs. ( 120]) and (|27|) plus the corrections terms in Eq. ( 1A-6|) lead to 



/ 



M 



I (NO) 



A 



yiv + y[v'ri y2v'-y' 2 v"r l (y' 3 v - y\V^A r^S 



-y& 



J 



M' D 



(IO) 



Vu 

A 



2/i u - V3V ^ J r t y 2 v + y 2 v n 

11, 11 10 \ \ 1 

y [Viv +ysv n y 2 v y 3 v 

V\V ~V3V -^j n y 2 v - y 2 v n 
j/5V - 2/3^"^) n 2/2^' + ZyWn 



-yzv- 



(A-14) 



to first order in r 1; in the NO and IO, respectively. As shown explicitly in the main text, 
the dynamics of the right-handed sector are relatively unaffected: the new entries in the 
first column of the Dirac mass matrices in Eq. ( 1A-14I) will induce mixing between the sterile 
neutrino and the fi and r flavors, but of the same magnitude as the original 9 e i, so that 
Q\ will not increase by that much [cf. Eqs. ( 135]) and ( 136]) ]. Thus the entire first column of 
M' D , suppressed by the mass M[ = O(keV), can be decoupled from the seesaw (assuming 
that \w[\ ^ \wi\)- In addition, corrections to U e5 in Eqs. ( 122]) and (128]) will also be small 
(see Sects. IIII C II and IIII C 2l for a discussion of those effects). 

The full 5x5 NLO neutrino mass matrix M' u 5x5 can now be constructed from the second 
and third columns of M' D and diag(M2,Mg), as in Eq. (119j) . Since we consider scenarios 
where NLO seesaw terms are important, we once again perform the full 5x5 diagonalization 
[cf. Eqs. (122]) and (128]) ]. including the new terms from higher-order operators in Eq. (1 A- 14j) . 



The matrix diagonalizing M' u 5x5 is explicitly given by 



| ^ TBM ° 3X2 | 

02x3 12x2 



(A-15) 
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where, to first order in n and second order in q, 

-^(x-PB + ^f) (1-P2)ei (P3-X)e 2 > \ 
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^[ 2e i + <(l + J R) e f] 
-^[ 2e |-^(l + i?) e ?] 











o \ 





3 £ 2 l-iV, , 
~2 e l 2°+ e l e 2 



-I / 



(A-16) 



in the normal ordering, where only first order terms in R ~ (9(10 x ) are kept [see Eq. 
and <r± = x ± p 2 — p 3 . The new mass eigenvalues are 

m\ = , 

m' 2 - |l - 3e? - f - i [9pl - 4p 2 ( x - p 3 ) - ( x - pa) 2 ] 6? 

-^[p2(l-36 2 )-< e |]} , 



m' 3 * {l - 2el + ( x - p 3 ) 2 (l - 3e|) - &±L{1 + 2R)e\ 

+ \[pI + 4 Mx ~ Ps) + 3(x - Ps) 2 ] - 2e|)j 

2, 3«)2 m (0) 



(A-17) 



m' 4 ~ W2u'X 2F2 — m 2 



l-3e 2 + ^(l- 



??? 



(o) 



m 5 ~ W3U 



{l-2 £ l + i(x " Pa) 2 - J [8(x - P3) 2 + (^) 2 R] 4} 



which corresponds to Eq. (124)) in the limit (XiPZiPz) ~~ 0- Here one can explicitly see that 
NLO seesaw corrections are controlled by whereas corrections from higher-order operators 
are controlled by \i P2 and p%. In those scenarios where the 6{ are negligible, i.e. scenario III, 
one could still have corrections from the latter. Those turn out to be small in the normal 
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ordering. 

In the inverted ordering, we have 



( ^(3p 2 + 4G) -^ M (3p 2 + a I + G) (l-p 2 )6i (2- X + p 3 )e 2 \ 
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(3 - 2 X + 2p 3 ) < 



V ^ 

to first order in x and second order in e,, where cr+ and G = (9(10 2 ) are defined in Eq. (1391) . 
In this case we cannot expand in G, in contrast to the NO case, where we expanded to first 
order in R. The new mass eigenvalues are 

m\ ~ mS 0) |l-6e|-i( X - p 3 ) (6 - <ri) + [8( X - Ps) + 1p\ + 4p 2 ( X - p 3 ) - 4(* - p 3 ) 2 ] e 

+^ N " (X " Ps) 2 ] G + [3p 2 + 4p 2 ( X - p 3 ) + (x ~ Ps) 2 } Gel + ^G 2 (l - 3ej) } , 
,^ ( _ _ _ 1 _ 



m' 2 ~ mf |l - 3(1 + 6p 2 )e 2 + 4p 2 + -L [27p 2 + 12p 2 ( X - p 3 ) + 
[6(1 + 2G)e 2 - G 2 (l - 6e 



( X -p 3 ) 2 ] G(l-3e 2 ; 
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In this case the corrections very much depend on the scenario concerned, since the value 
of the €i terms can give cancellations. However, the correction to \U e 2\ 2 constrains the 
parameters \i Pi an d p3 to be small (see discussion in the main text), and since G always 
occurs together with one of the three parameters the effect of G = O(10 2 ) will always be 
suppressed. In the end we always have enough parameters to fit the mass eigenvalues to the 
data. 
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